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Abstract 

A search for excited state transitions of double beta decay in 110 Pd and double electron capture in 
102 Pd has been performed in the HADES underground laboratory with two low background HPGe 
detectors in sandwich configuration. After an exposure of 35.92 kg-d, no signal was found. The 
frequentist spectral analysis resulted in lower half-life limits of 1.98 • 10 20 yr and 1.72 • 10 20 yr (95 % 
CL) for the first + and 2 + excited state in 110 Pd respectively. This is an improvement by more 
than a factor of 3 with respect to previous measurements. In 102 Pd, the lower half-life limit could 
be improved to 5.95 • 10 18 yr (95 % CL) for the first + excited state. Furthermore, first experi- 
mental lower half-live limits are found for all possible excited states in the 110 Pd and 102 Pd systems. 



Structured Abstract: 

Background: Excited state transitions in double beta decays are a powerful tool to validate and 
tune calculations of nuclear matrix elements. 

Purpose: The experimental lower half-life limits for double beta decays of 110 Pd and 102 Pd into 
the excited 2f and 0^ are improved. Furthermore, first limits of transitions into the 2^, 0^ and 
2g~ states are published for 110 Pd as well as a first limit for the 2^~ state transition in 102 Pd. 

Methods: The Pd sample was measured with two HPGe detector in sandwich configuration in 
the HADES underground laboratory during 44.77 d of life time. The analysis is performed with 
the frequentist Fieldman Cousisns method. 

Results: Lower half-life limits of 1.98 • 10 20 yr and 1.72 • 10 20 yr (95% CL) have been found for 
the first + and 2 + excited state in 110 Pd respectively. This is an improvement by more than 
a factor of 3 with respect to previous measurements. In 102 Pd, the lower half-life limit could 
be improved to 5.95 • 10 18 yr (95 % CL) for first + excited state. Furthermore, first experi- 
mental lower half-live limits are found for all possible excited states in the 110 Pd and 102 Pd systems 

Conclusions: Previous half-life limits have been improved and experimental results were obtained 
for all theoretical calculations of palladium double beta decays into excited states. 
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I. INTRODUCTION 



The field of neutrino physics has made groundbreaking progress in recent years. Oscilla- 
tion experiments studying neutrinos coming: from the sun l|-|3|, the atmosphere 4j], nuclear 



reactors 



-[8] and in accelerator beams |9Nll| have found compelling evidence for flavor os- 
cillation in the lepton sector. This changes the long believed assumption of neutrinos being 
mass-less particles. The implication from neutrino oscillation, that at least two neutrino 
mass eigenstates have a non-zero rest mass, does not allow fixing an absolute mass-scale and 
leaves two mass hierarchy scenarios open. Advanced oscillation experiments try to identify 
the hierarchy scenario by using oscillation effects in the earth or in dense stellar matter 12] 
3ut are not able to determine the absolute mass scale. This has to be done via beta-decay 
13, Q, cosmology [lj| or neutrinoless double-beta decay: 



(Z,A) -> (Z + 2,A) + 2e~ (0uf3f3) , (1) 

which would have a high sensitivity to determine the mass scale of neutrinos. This process 
violates total lepton number by two units and thus is not allowed in the Standard Model. 
Furthermore, it is the gold plated process to distinguish whether neutrinos are Majorana or 
Dirac particles. 

For the 0uj3j3 process to exist, it is necessary to match the helicities of the intermediate 
neutrino states which is easiest done by introducing a neutrino mass. This mass is connected 
to the experimentally observable half-life via 

(T^y 1 = G°»(Q, Z) ■ \M°J T - M°f ■ (^A\ 2 , ( 2 ) 

\ m e J 

where (jn Ue ) is the effective Majorana neutrino mass given by the coherent sum over the 
virtual electron neutrino mass eigenstates (m Ve ) = \ U^m^ with U ei as the lepton flavor 
mixing matrix, G 0u (Q,Z) is a phase space factor and Mq u t — Mp u describes the nuclear 
transition matrix element. The experimental signature is the emission of two electrons with 
a sum energy corresponding to the Q-value of the nuclear transition. A potential evidence 
ias been claimed for the 0u(3f3 mode of 76 Ge with T°Y 2 = 2.23l{j;^ x 10 25 yr at 90% CL 
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Experimentally observed in 11 nuclides 18|,ll9| is the SM process of neutrino accompanied 
double beta decay: 



[Z, A) -> (Z + 2, A) + 2e~ + 2u e (2u/3/3) 



(3) 



which is expected with half-lives around 10 20 yr depending on the Q-value. For recent 
reviews see 20|]. An alternative process is the double positron decay in combination with 
electron capture (EC). Three different decay modes can be considered: 



(Z,A)^(Z-2,A) + 2e + + (2u e ) (f3 + /3+ ) (4) 

e- + (Z,A) ^(Z - 2, A) + e + + (2u e ) (/? + /EC) (5) 

2e~ + (Z,A) ->(Z-2,A) + (2u e ) (EC/EC) (6) 

Decay modes containing a positron have a reduced Q-value as each generated positron 
accounts for a reduction of 2 m e c 2 and thus can only occur in nuclides with sufficient energy 
difference to the daughter nuclide. The full energy is only available in the EC/EC mode 
and makes it the most probable one; however, it is also the most difficult to detect, only 
producing X-rays instead of 511 keV gammas. 



Another set of searches focuses on excited state transitions in double beta decays. This is 
experimentally interesting because the event topology is enhanced by de-excitation gammas 
which are accessible with gamma ray spectroscopy. However, without information about the 
other final state particles, this technique cannot distinguish between the Ovfifi and the 2v/3/3 
mode; hence, the deduced half-lives are valid for both. The investigation of 2v/3/3 modes 
into excited states provides information on nuclear structure which is valuable for matrix 
element calculations. The understanding and tuning of parameters in matrix elements for 
2v/3{3 modes is imperative for translating the measured half-life of Ou/3/3 experiments into 
a Majorana neutrino mass (Eg. [2D- So far only transitions to the first excited + state 
have been observed in 100 Mo 2l| and 150 Nd 22|. The searches described in this paper are 
searches for excited state transitions in palladium with gamma ray spectroscopy. 
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II. DOUBLE BETA DECAYS IN PALLADIUM 



The element under study is palladium with the isotopes of interest 110 Pd and 102 Pd. 
Among the 35 isotopes expected to undergo f3~f3~ decay, 110 Pd has the second highes natu- 
ral abundance with 11.72%. Recently, the Q-value was remeasured to 2017.85(64) keV 23] 
and places 110 Pd among the 11 (3~f3~~ isotopes with a Q-value larger than 2000 keV. Two 



measurements of 110 Pd have been performed in the past in 1954 24j and more recently in 
2011 25j. The latter measurement was the first to investigate excited states in palladium 
and is the direct predecessor of this search. There exist many theoretical calculations for 
transitions into the ground state and into the 2±, 0^, 2 J, 0^ and 2% excited states to 
which the experimental limits can be compared. The existing experimental and theoretical 
half-life limits ares summarized in Tab. [I] This paper has the aim to provide experimental 
information about every decay mode into excited states in 110 Pd and 102 Pd that has been 
investigated theoretically. Furthermore, 110 Pd is an excellent candidate to probe the single- 
state dominance hypothesis for 2^/3/3-decay, i.e. that only the lowest lying intermediate 
l + -state will contribute to the nuclear transition matrix element describing its 2i//3/3-decay. 



The second isotope 102 Pd has a Q-value of 1172keV, a natural abundance of 1.02 % and 
is able to decay via EC/EC and /3 + /EC. The only experimental half-life limit was quoted 
in {25I and no theoretical calculation have been published up to date. A summary can be 
found in Tab. HT1 



The search is based on gamma spectroscopy, hence only gamma lines are considered in the 
event topology. Each excited state transition is followed by an unique set of decay branches 
and gamma cascades which are illustrated in Fig. [1] and [2] 



III. EXPERIMENT 



A. Measuring Setup 

The measurements were performed in the High Activity Disposal Experimental Site, 
HADES, underground laboratory on the premises of the Belgian Nuclear Research Centre 
SCK-CEN in Mol. The underground laboratory is located at a depth of 225 m inside a 
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TABLE I. Experimental and theoretical half-life limits for various /3/3 decay modes in Pd. The 
columns show from left to right the theoretical model the quoted half-life, the reference and the 
year of publication. Abbreviations denote: PHFM - Projected Hartree-Fock-Bogoliubov, SSDH - 
Single-State-Dominance Hypothesis, SRPA - Second quasi Random Phase Approximation, OEM 
- Operator Expansion Method, QRPA - Quasi Random Phase Approximation, SSD - Single State 



Dominance and pnQRPA - proton-neutron quasiparticle random-phase approximation. 



Exp / Th model 


Lower Limit Ty 2 [yr] 


Reference 


Year of Publication 


110 Pd Ground State Transition 


Exp 


1 ■ 10 17 (68 % CL) 


[24] 


1952 


PHFM 1.41 ■ 10 20 and 3.44 ■ 10 20 ^ 


[26] 


2005 


SSDH 


1.75 ■ 10 20 


[27] 


2000 


SSDH 


1.2 - 1.8 ■ lO 20 ^ 


[28] 


1998 


SRPA 


1.6 ■ 10 20 


[29] 


1994 


OEM 


1.24 ■ 10 21 


[30] 


1994 


QRPA 


1.16 ■ 10 19 


131] 


1990 


SSD 


1.2 ■ 10 20 


32] 


2005 


pnQRPA 1.1 ■ 10 20 and 0.91 ■ 10 20 ^ 


Tool 

33] 


oni 1 
ZU1 1 


110 Pd 2+ Excited State Transition 


@ 657.76 keV 






Exp 


4.40 ■ 10 19 (95 % CL) 


[25] 


2011 


SSD 


4.4 ■ 10 25 


[32] 


2005 


SRPA 


8.37 • 10 25 


[29] 


1994 


pnQRPA 


1.48 ■ 10 25 


[31] 


2007 


pnQRPA 0.62 ■ 10 25 and 1.3 ■ 10 25 ^ 


[33] 


2011 


110 Pd 0^ Excited State Transition 


@ 1473.12 keV 






Exp 


5.89- 10 19 (95% CL) 


[25] 


2011 


SSD 


2.4 ■ 10 26 


[32] 


2005 


pnQRPA 


4.2 ■ 10 23 and 9.1 • 10 23 ^ 


[33] 


2011 


110 Pd 2+ Excited State Transition 


@ 1475.80 keV 






SSD 


3.8 ■ 10 31 


[32] 


2005 


pnQRPA 


11 ■ 10 30 and 7.4 ■ 10 30 ^ 


[33] 


2011 


110 Pd 0^ Excited State Transition 


@ 1731.33 keV 






SSD 


5.3 • 10 29 


[32] 


2005 


110 Pd 2% Excited State Transition 


@ 1783.48 keV 






SSD 


1.3 ■ 10 35 


[32] 


2005 



a gA = 1-25 and 1.0 respectively 

b Different experimental input for calculations 

c For Woods-Saxon Potential and adjusted base respectively (see 



33| for details) 
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TABLE II. Experimental half-life limits for various ECEC and /3 + EC decay modes in Pd. 



Exp / Th model Lower Limit Ti/ 2 [yr] 


Reference 


Year of Publication 


102 Pd 2f Excited State Transition @ 475.10kcV 


Exp 2.68- 10 18 (95% CL) 


[25] 


2011 


102 Pd 0+ Excited State Transition @ 943.69 kcV 


Exp 7.64- 10 18 (95% CL) 


[25] 


2011 



24.6s 



110 



Pd 



no 



Rh 



Qbb=2018keV 



\0+2 






\2+ 2 








J* 


\0+i 
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■OkeV 



FIG. 1. Pd level scheme of investigated decay modes. Nuclear data from 



Boom clay formation and has a flat overburden that amounts to roughly 500 m w.e. 36 j 



The detector setup consists of two HPGe detectors in a sandwich configuration with 



integrated muon veto panels on the top [37|. It is shown in Fig. [31 The sample is placed 
between the top (Ge-7) and bottom (Ge-6) detector. The distance between the two can 
be adjusted in order to maximize the solid angle acceptance and the detection efficiency. 
Ge-6 is a p-type HPGe detector with 80% efficiency and 0.9 mm dead layer in a 102 mm 
cryostat with Cu endcap whereas Ge-7 is an extended range p-type HPGe detector with 
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FIG. 2. 102 Pd level scheme of investigated decay modes. Nuclear data from |35l ]. 



90 % efficiency and 0.3 /im dead layer in an Al cryostat. The characteristics of Ge7 makes 
it suitable for the detection of low energetic X-rays while the configuration of Ge6 has the 
advantage of reducing the background and X-ray coincidences. The shielding consists of 
an outer layer of 14.5 cm 20Bq/kg ( 210 Pb) lead, an intermediate layer of 4.0 cm 2.4Bq/kg 

3.5 cm electrolytic copper with less than 



( 210 Pb) low activity lead and an inner layer o 
15 /iBq/kg 60 Co and less than 20 /iBq/kg 228 Th 



37]. 



The data acquisition is two-fold. The main DAQ is the IRMM self-fabricated DAQ2000 
multi-parameter system that is recording events in list mode from the two HPGe detectors as 



we 



1 as from the two muon panels and enables coincidence analysis in the ROOT framework 
]. Additionally, a standard GENIE DAQ system was used in histogram mode for each 
HPGe detector redundantly. The DAQ2000 was only operative for a reduced measuring 
time with limited sample exposure; thus the analysis in this work is performed with the 
GENIE DAQ and without muon veto or detector coincidence. 



The total background rate in the germanium detectors was previously measured with 
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HI electrolytic copper 

ultra low background lead 
low background lead 



FIG. 3. IRMM germanium detector sandwich setup. 



992cts/d in an energy range of 40..2400keV of which 124cts/d were identified as muon 



events 



37]. 



B. Palladium Sample 



802.35 g of irregular shaped 1 mm x 1 cm 2 plates of palladium were placed inside a measur- 
ing container of 70 mm diameter and 50 mm height. The plates are piled inside the container 
as dense as possible and an effective density of 10.2 g/cm 3 was calculated. The palladium is 
approximated with a homogeneous distribution and the effective density in the simulations 
for determining the detection efficiency Recently, the sample was purified by C. HAFNER 
GmbH + Co. KG in 2010 to a certified purity of > 99.95 % which lowered the continuous 



background in the peak regions by approximately 20 % [25j. In order to avoid radioisotopes 
produced by cosmic ray spallation, the palladium was kept underground and exposed only 
18 d during purification in 2010 and 3 d for transport in fall 2011. A picture of the palladium 
sample before and after purification is shown in Fig. HI 
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FIG. 4. Palladium sample before (left) and after (right) purification inside the measuring con- 
tainer. 

IV. ANALYSIS 

A. Stability Check 

The histogram data of the GENIE DAQ is separated into individual runs of roughly 24 h 
for each of the two detectors Ge6 and Ge7. This enables the stability check of the DAQ 
system over the extended period of measurement and the removal of individual runs. 

The stability was checked by plotting the count rate and peak centroid of the background 
peaks at 609 keV and 2614 keV as a function of time (spectrum number). Although counting 
statistics prevented detailed analysis of the short term stability, it was clear that there was 
no measurable energy drift during the measurement period. In addition, quality controls 
with a point source containing 60 Co, 137 Cs and 241 Am were carried out before, after and once 
in-between the measurement. The average 222 Rn activity concentration in the laboratory 
during the data taking was measured to 9.5Bq/m 3 and at no point higher than 25Bq/m 3 . 
No correlation to the background count rate was observed mainly due to effective Rn-removal 
by minimizing empty space inside the shield and flushing with N 2 . The count rate remained 
stable over the 32 runs with one exception; in the first two runs an increase in total counts 
and counts from 208 T1 was observed, however not from 214 Bi. This behavior was cross checked 
with additional visible 7-lines from 583.19keV 208 T1 and 239.63 keV 212 Pb, and 1764.49 keV 
214 Bi representing the 232 Th and 238 TJ decay chains respectively. The increased count rate 
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of the 232 Th chain in the first days is only seen in the lower Ge6 with a Cu endcap and not 
in the upper Ge7 with an Al endcap. One possible explanation is that the short half-life of 
220 Rn (55.6 s) implies that it is not flushed out by boil-off nitrogen and that its daughters 
preferentially stick to the copper surface of the lower detector (Ge6) than to the aluminum 
surface of the upper detector (Ge7). 222 Rn with longer half-life (3.8 d) will be flushed out of 
the shield before a significant amount of daughters are produced. 

Consequently, the first two runs are removed from the analysis resulting in a total of 30 
runs with 44.77 d of good data out of 32 runs with 46.49 d total data which translates into 
a total exposure of 35.92 kg-d. 

B. Data Processing 

For the final analysis a single energy spectrum is used for which all individual runs are 
summed: In a first step, all selected runs of one detector are combined with the same 
energy calibration. In a second step, the single detector spectra are rebinned into a common 
binning of 0.5keV/bin and a common energy range from 20..2720keV. Additionally, the 
Ge6 spectrum is scaled to the life-time of the Ge7 spectrum which becomes the common 
life-time of the sum spectrum. The commonly binned and scaled spectra are added. These 
steps result in a non-integer sum spectrum that denotes the count per bin in the life-time of 
Ge7. The difference in life-time between the two detectors is less than 0.5 % for the selected 
runs. The combined spectrum of Ge6 and Ge7 together with a background spectrum of 
13.62 d is shown in Fig. |5j 

C. Background Investigation 

The radiopurity of the palladium sample was assessed during the 7-spectrometry mea- 
surement performed in HADES. In the measured spectrum, the major 7-lines emitted by 
natural radionuclides belonging to the 238 U, 232 Th chains and 40 K as well as the 60 Co lines 
are visible. The palladium spectrum is compared to the background spectrum measured 
without a sample. The background peak count rate is subtracted from that of the palla- 
dium sample. The result, if positive, is then used for the evaluation of the activity due to 
the impurities in the sample. In case of a negative result, a decision threshold is calculated 
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FIG. 5. Sum spectrum of Ge6 and Ge7 for selected data set (44.77 d) in red and background 
spectrum (13.62 d) in gray. The spectra are shown in common 10 keV bins between 20 and 2720 keV. 



according to 39J. The results are reported in Tab. IIIIl Decision thresholds are also cal- 
culated for the following radionuclides: 102 Rh (Ti /2 = 207.3d), 102m Rh (Ty a = 3.742 yr) and 
110m Ag (Tu 2 — 249.76 d). The reason is the possible interference with the search for 110 Pd 
and 102 Pd isotopes, because of the emission of 7-lines from the same excited daughter states. 
No presence of these radionuclides is found, as reported in Tab. IIVI 



D. Peak Finding 

The analysis is a peak search on the detector sum spectrum and either retrieves the 
number of counts in a respective peak or states an upper limit of counts according to a level 
of confidence. 

In an experiment with non negligible background, the background can fluctuate upwards 
or downwards. The sensitivity of an experiment is then defined as a resulting signal which 
originates from a 1 a upward fluctuation of the background. This can be calculated before 
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TABLE III. Massic activities (in mBq/kg) of radioimpurities detected in the Pd sample. 









Decision Threshold 


Weighted Mean 


Nuclide 


E 


Massic Activity 


(a = 95 %) 


Massic Activity 




[keV] 


[mBq/kg] 


[mBq/kg] 


[mBq/kg] 


214 Pb 


295.22 


1.9 ± 1.0 


1.4 


1.4 ± 0.4 




351.93 


1.3 ± 0.5 


0.6 




214 Bi 


609.32 


1.9 ± 0.4 


0.4 


1.9 ± 0.4 




1120.29 


2.0± 0.8 


0.9 






1238.11 


— 


2.2 






1377.67 


— 


2.7 






1764.54 


— 


3.2 




21 °pb 


46.54 




414.3 




228 Ac 


911.20 


— 


0.5 






968.97 


— 


0.9 




212 Pb 


238.63 


— 


0.7 




208 rp^ 


583.19 


— 


0.6 






2614.51 




0.3 




40 K 


1460.82 




1.0 




137 Cs 


661.66 




0.2 




60 Co 


1173.23 




0.2 






1332.49 




0.1 





performing the experiment if the background is known. In case of the observation of a 
downward fluctuation, which formally results in negative signal counts, the signal is usually 
set to zero and the sensitivity is quoted as an upper limit of the counts. On the other 
hand, when using classical uncertainties on the observed downward fluctuated counts, it is 
possible that in some cases even the upper limit its negative. In these cases it results in a 
poor coverage of the quoted confidence level at best and in an unphysical negative result at 
worst. Additionally there is a discontinuity in the coverage when crossing from a two-sided 
confidence interval definition to a one-sided one. 
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TABLE IV. Decision thresholds for direct 7-background of DBD intermediate nuclei. 



Decision Threshold 



Nuclide E {a = 95 %) 

[keV] [mBq/kg] 

110m Ag 1384.30 0.5 

1505.04 1.0 

102 Rh 475.05 0.4 

102m Rh 631.28 0.3 

697.49 0.3 



All these problems are addressed by the method of Feldman and Cousins in their paper 



40| . They use a Neyman construction of a confident belt and an ordering principle based 



on likelihood ratios. The advantages are a physical yield in all background situations i.e. 
positive upper count limit and avoiding discontinuities in the coverage while crossing the 
statistical interpretation from a non-observation to an observation, i.e. from a one-sided to 
a two-sided confidence interval. The confidence intervals of the Feldman Cousins method 



are believed to have a better coverage for small numbers than gaussian ones 40]. 



No prominent peak structures are observed in the signal region and the results of the 
analysis are upper limits only for the peak counts. In order to obtain a numerical value, all 
bins within a peak are combined into a single analysis bin that covers at least the FWHM 
of the peak. The real signal fraction coverage is calculated as the gaussian peak area in the 
analysis bin and depends on the actual binning of the spectrum. The peak background is 
fitted with a constant functions defined ±30 keV around the peak energy excluding a window 
of ±5 keV around the peak. In case of prominent background peaks in the side bands, they 
are included in the background function as gaussians. This was done for the background 
peaks at 609.31 keV, 1120.29 keV and 1764.49 keV from 214 Bi, at 238.63 keV from 212 Pb, at 
1461.83 keV from 40 K and at 661.66 keV from 137 Cs. In case of the 661.66 keV peak which is 
close to the 657.76 keV peak in the 110 Pd system, the flat background function was defined 
closer than ±5 keV into the signal region to improve the background estimation. 
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The observed counts in the analysis bin are compared to the expected background with 
the ROOT class TFeldmanCousins, which returns the lower and upper bound of the signal 
confidence interval according to a specified confidence level that is set to 95 % in this work. 
All investigated peaks show a lower count limit of zero; this is in agreement with a non- 
observation of the peak. To account for the incomplete coverage of the peak area by the 
analysis bin, the upper count limit is divided by the fraction of coverage and thus adjusts 
the upper count limit in a conservative way. 



An illustration of the technique is shown in Fig. to with the energy spectrum in the 
solid black line, the background function in the dashed blue line, the peak fraction marked 
as the solid red area and the gaussian signal peak as it appears with the Feldman Cousins 
upper limit as the red solid line. 



39] 



The Feldman and Cousins limits are cross checked with the ISO standard methods 
and agree better than within a factor of two with each other; this can be explained by the 
different treatment of statistical background fluctuations. 



E. Monte Carlo Simulation for 7-line Efficiencies 

The full energy peak efficiencies are determined using Monte Carlo simulations with 
the EGS4 software. The models of the detectors were first determined from manufacturer 
data and using information from radiography. Thereafter the dead layer thicknesses were 
adjusted in the model to agree with measured FEP efficiencies from point sources within 
3 %. The final model has been validated using volume sources of similar size as the Pd 
source in this study. Each decay branch is simulated separately with information from 3s| 
and the calculations involved all the cascading 7-rays of each branch so that the resulting 
FEP efficiency is inherently corrected for the coincidence summing effect. X-ray coincidences 
and the angular correlations are neglected in the simulations and it was assumed that the 
activity was homogeneously distributed in the whole volume of the sample. 
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V. RESULTS 



All 7-lines participating in a 7-cascade are investigated and a half-life was calculated for 
each. Intrinsically, the calculation of limits is influenced by statistical fluctuations in the 
experimental spectrum, hence the largest calculated limit for an excited state transition is 
quoted as the half-life limit of this transition. A summary of all investigated 7-lines can 
be found in Tab. |Vj Quoted are the emission probability, the detection efficiency including 
summation effects, the upper count limit in the spectrum and the calculated half-life. The 
selected half-life for each transition are summarized in Tab. I VI I The peak regions for the 
7-lines originating from the favored Of transitions are shown in Fig. and [7] for 110 Pd and 
in Fig. M and M for 102 Pd. 




FIG. 6. The peak region around the 657.76 keV 7-line in 110 Pd and the 661.66 keV 137 Cs 7-line. 
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TABLE V. Experimental results for each decay mode and 7-line. The columns from left to 
right denote the decay mode, the 7-line energy, the emission probability in % and the detection 
efficiency in %. The last two columns show the upper signal count limit and the deduced lower 
half-life limit at 95 % CL. 



Decay Mode 


■y-line Energy [keV] 


Emission Probability 


Detection Efficiency 


Signal Count Limit 


1 1/2 -Limit [yrj 


110 "D J 0+ IP 

r& OOl.lOKCV 


pki ip ~\ r ^\r 

OO ( . ( KC V 


1UU /o 


4. / U /o 


1 O A 


1 70 1 n20 
1. ( Z • 1U 


110 pj n+ 1,170 1 9 V pV 


Q1 c: OQ \ra\7 

010.00 Ke v 


1 no % 


o.o*± /o 


Si l 
0.1 


1 qq . 1 n20 




657.76 keV 


100% 


3.94% 


12.4 


I.44 . 10 20 


110 Pd 2% 1475.80 keV 


1475.80 keV 


35.25 % 


1.32% 


11.5 


5.17 ■ 10 19 




818.02 keV 


64.75 % 


2.40% 


16.3 


6.67 ■ 10 19 




657.76 keV 


64.75 % 


2.53% 


12.4 


9.26 ■ 10 19 


110 Pd 0% 1731.33 keV 


1073.7 keV 


86.73 % 


1.89% 


10.1 


8.50 ■ 10 19 




00 / . 1 keV_ 


OK Q O 0/ 

yo.32 70 


3. to 70 


lz.4 


1.38 ■ 10 




255.49 keV 


13.27% 


0.36% 


25.3 


6.46 ■ 10 18 




1475.80 keV 


4.68 % 


0.12% 


11.5 


4.87 ■ 10 18 




818.02 keV 


8.59 % 


0.24% 


16.3 


6.63 ■ 10 18 


110 Pd 2+ 1783.48 keV 


1783.48 keV 


21.57% 


0.88 % 


6.2 


6.45 ■ 10 19 




1125.71 keV 


78.43 % 


2.48% 


12.0 


9.41 ■ 10 19 




657.76 keV 


78.43 % 


2.99% 


12.4 


1.09 ■ 10 20 


102 Pd 2+ 475.10 keV 


475.10 keV 


100% 


5.09% 


33.7 


5.95 ■ 10 18 


102 Pd 0+ 943.69 keV 


468.64 keV 


100% 


4.32 % 


29.3 


5.81 ■ 10 18 




475.10 keV 


100% 


4.31 % 


33.7 


5.04 ■ 10 18 


102 Pd 2+ 1103.05 keV 


1103.05 keV 


37.11% 


1.60% 


13.5 


4.66 ■ 10 18 




627.94 keV 


62.90 % 


2.54% 


11.7 


8.55 ■ 10 18 




475.10 keV 


62.90 % 


2.67% 


33.7 


3.13 ■ 10 18 



a This 7-line is part of two sub-branches starting from the same excited state. 
CONCLUSION 

A palladium sample has been investigated for double beta decay transitions into excited 
states in the low background laboratory HADES. Lower half-life limits could be improved 
for the 0^ and 2^ transitions in 110 Pd and 102 Pd and first limits are established for all 
possible higher energetic excited state transitions. The best limit could be set for the 110 Pd 

transition with a half-life larger than 1.98 • 10 20 yr. The largest improvement compared 
to previous results is archived for the 110 Pd 2\ transition with 1.72 • 10 20 yr which is an 
improvement by a factor of 3.9. For the 102 Pd system, the improvements are smaller due to 
upward fluctuations of the background in the peak region of the 475.10 keV and 468.64 keV 
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TABLE VI. Summary of measured half-life limits for all Pd and Pd double beta decay 
excited state transitions. 



Decay Mode 



Ty 2 Limit [yr] (95 %] 



110 Pd 2+ 657.76 keV 
110 Pd 0+ 1473.12 keV 
110 Pd 2+ 1475.80 keV 
110 Pd 0^ 1731.33 keV 
110 Pd 2+ 1783.48 keV 



1.72 • 10 20 
1.98 • 10 20 
9.26 • 10 19 
1.38 • 10 20 
1.09 • 10 20 



102 Pd 2+ 475.10 keV 
102 Pd 0+ 943.69 keV 
102 Pd 2+ 1103.05 keV 



5.95 • 10 18 
5.81 • 10 18 
8.55 • 10 18 



I 



12 



10 



'J] 



Bg fit: 25.2 cts/ROI (5.0 cts/bin) 
Obs cts: 21.5 cts/ROI (4.3 cts/bin) 
FC upper limit: 8.4 cts (95% CL) 
ROI peak fraction: 0.85 
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FIG. 7. The peak region around the 815.33 keV 7-line from 110 Pd 0?". 



7-lines. 



Possible improvements of the search of double beta decays in palladium is the consider- 
ation of X-rays. This requires a different geometric assembly of the palladium plates e.g. in 
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FIG. 8. The peak region around the 475.10 keV 7-line from Pd 2J~. 




FIG. 9. The peak region around the 468.64 keV 7-line from Pd 0^ . 

a layer around an n-type HPGe detector with thin dead layer. Further improvement can be 
achieved considering 7— coincidences using the multi-parameter DAQ system. This would 
also reduce the muonic background. The intrinsic massive activity of the palladium sample 
was determined to be 1.7mBq/kg for the 238 U chain and below the detection threshold for 
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the 232 Th chain. Further purifications are not expected to yield significant improvement. 
However, accumulating storage underground will reduce the general background from cosmic 
activated radio isotopes in the palladium sample and the measuring system. 
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